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The VERITAS array of 12-m atmospheric-Cherenkov telescopes in southern Arizona is one of 
the world's most-sensitive detectors of very-high-energy (VHE, E >100 GeV) gamma rays. More 
than 50 extragalactic sources are known to emit VHE photons; these include blazars, radio galaxies, 
and starburst galaxies. Blazar observations are one of the VERITAS Collaboration's Key Science 
Projects. More than 400 hours per year are devoted to this program and ~ 100 blazars have already 
been observed with the array, in most cases with the deepest ever VHE exposure. These obser- 
vations have resulted in 21 detections, including 10 VHE discoveries, all of them with supporting 
multiwavelength observations. Recent highlights from VERITAS extragalactic observation program 
and the collaboration's long-term blazar observation strategy are presented. 
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I. THE POPULATION OF VHE BLAZARS 

Blazars are active galactic nuclei (AGN) hosting a 
supermassive black hole, with relativistic jets point- 
ing close to the line of sight to the observer. The 
small viewing angle of the jet makes it possible to ob- 
serve strong relativistic effects, such as a boosting of 
the emitted power and shortening of the character- 
istic time scales. Blazars are typically identified in 
radio and optical surveys, their defining characteris- 
tics being radio loudness, flat spectrum and compact 
morphology in the radio band, some degree of optical 
polarization, and fast, large-amplitude variability. 

Since the first detection of a blazar in the very- 
high-energy gamma-ray band (VHE, E > 100 GeV) in 
1992 by the Whipple Collaboration [l|, ground-based 
imaging Cherenkov telescopes have discovered up to 
45 VHE-emitting blazars This number is growing 
fast, with half of these detections reported in the last 
three years. VHE blazars populate the whole sky. Fig- 
ure[T]shows the distribution of VHE blazars in galactic 
latitude and declination. The only noticeable asym- 
metry is that two thirds of the sources are located in 
the northern half of the sky. In the redshift distribu- 
tion of VHE blazars (see Figure [2]) most objects occur 
at z < 0.2, with the most distant known VHE source 
being 3C 279 at z = 0.536 Q. The pair-production of 
VHE photons with UV to infrared photons from the 
extragalactic background light (EBL) limits the trans- 
parency window for gamma-rays at large redshifts. Fi- 
nally, a significant fraction of VHE blazars (10 out of 
45) have unknown redshift, or published redshifts that 
have been reported as uncertain. 

Blazars display a characteristic two-component 
spectral energy distribution (SED, see Figure[3]). The 
low-frequency continuum is usually attributed to syn- 
chrotron emission by relativistic electrons inside the 
jet magnetic field. VHE observations of blazars sam- 
ple the high-frequency peak of the SED. Most theo- 
retical models assume that the gamma-ray emission 



comes from inverse-Compton (IC) scattering of low- 
energy photons off of a population of relativistic elec- 
trons. The source of photons for the IC process can be 
either synchrotron photons emitted by the same elec- 
tron population [so-called synchrotron self-Compton 
models (SSC), e.g., @], or a population of photons ex- 
ternal to the jet [external Compton (EC) models, e.g., 
Q . Another family of models includes emission from 
hadrons in the jet [e.g., establishing AGN jets as 
candidates for the origin of the highest energy cosmic 
rays [see [Toj . 

Table Q] shows the average properties of VHE- 
detected blazars divided by blazar class. Most de- 
tected sources are of the BL Lac class, with the 
only exception of three flat-spectrum radio quasars 
(FSRQ). VHE-emitting FSRQs exhibit very soft spec- 
tral index and are only detected during very bright 
flaring episodes @, Gjf . FSRQs are also the most 
distant sources detected in the VHE range. High- 
frequency-peaked BL Lacs (HBL) are the dominant 
subclass of VHE-detected blazars, accounting for 
three quarters of the total source count. HBLs are 



Class # of sources {z) a (T) b 



t v a 



HBL 
IBL 
LBL 

FSRQ 
Total 



33 
5 
4 
3 

45 



0.12 
0.12 
0.14 
0.44 



3.2 
3.7 
3.2 
3.9 



< 5 min, [11, 12} 
~lday, [13] 
< 20 min, [14] 
~ 10 min, [15] 



TABLE I: Average properties of all VHE-detected blazars 
divided by class: high, intermediate, and low-frequency- 
peaced BL Lacs (HBL, IBL, LBL, respectively) and flat 
spectrum radio quasars (FSRQ). Data taken from 
a Average redshift, excluding the sources where the red- 
shift is unknown or uncertain. 
6 Average photon index in the VHE range. 
c Shortest observed variability time scale in the VHE band, 
with the corresponding reference. 
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FIG. 1: Top: All-sky map of all VHE-detected blazars in galactic coordinates from [jj. The shaded area indicates 
the region observable with VERITAS. Bottom: Distribution of VHE blazars in galactic latitude (left) and celestial 
declination (right). The corrsponding distributions are shown for candidate VHE blazars || and .Fermi-detected blazars 
[3] for comparison. 



defined by having the peak of the low-energy (or 
synchrotron) component (v syn c) a t frequencies higher 
than 10 15 Hz. Observations of hard-spectrum, mid- 
redshift HBLs have been used in recent years to place 
strong constraints on the density of the EBL [3 [I?} , 
to study the formation of pair cascades in the propaga- 
tion of TeV photons through the intergalactic medium 
and to derive constraints on the strength of intergalac- 
tic magnetic fields [18-20]. Intermediate-frequency- 
peaked BL Lacs (IBL, 10 14 Hz < v sync < 10 15 Hz) 
and low-frequency-peaked BL Lacs (LBL, v sync < 
10 14 Hz) are typically detected during high flux states, 
with low gamma-ray states falling at the level or be- 
low the sensitivity of current instruments. The SED 
of VHE-detected IBLs and LBLs often shows evidence 
for a significant contribution of EC processes to the 
total gamma-ray output [e.g., see Figure [3J. 

The fastest flux variability in blazars has been ob- 
served in the VHE range, with individuals of most 
blazar subclasses exhibiting sub-hour flux variations 
(see Table [TJ . The observed flux-doubling times indi- 



cate compact emission regions with size < 10 -3 pc 
[22| . Such compact regions suggest that the emis- 
sion could originate from the base of the jet, in con- 
tradiction with the downstream dissipation scenarios 
favored by simultaneous radio and optical observa- 
tions of blazars during gamma-ray flares 0, [24| . Al- 
ternatively, very compact emission regions embedded 
within the large-scale jet [25, 26], jet recollimation [27j 
or magnetic reconnection events [28j could accommo- 
date fast variability with far dissipation scenarios. 



II. VERITAS BLAZAR OBSERVATIONS 

VERITAS is an array of four imaging atmospheric 
Cherenkov telescopes located at the Fred Lawrence 
Whipple Observatory (FLWO) in southern Arizona 
(31°40'N, 110°57'W, 1.3km a.s.L). It combines a 
large effective area over a wide range of energies 
(100 GeV to 30 TeV) with an energy resolution of 15- 
25% and an angular resolution of less that 0.1°. The 
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FIG. 2: Redshift distribution of all VHE-detected blazars, 
including the distribution for each blazar class. Data taken 
from Q. Blazars with unknown redshift include objects 
without a quoted redshift in 0] and blazars with published 
redshift values that have subsequently been reported as 
uncertain. 
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FIG. 3: Broadband SED of 3C 66A during the 2008 Octo- 
ber multiwavelength campaign [g|. The dotted line shows 
an SSC+EC model assuming a redshift of z = 0.3. 



high sensitivity of VERITAS allows the detection of 
sources with a flux of 0.01 times that of the Crab 
Nebula in about 25 hours. The standard VERITAS 
analysis methods are described in (2t| [3(| ■ 

Table [TT] shows the list of 21 blazars detected in 
the VHE band with VERITAS since it started regu- 
lar observations in 2007. Ten of these detections are 
first-time discoveries in the VHE band. VERITAS 
typically accumulates ~ llOOh/year of observations, 
with blazar observations adding up to ~ 400h/year 
[3ll ] . The VERITAS blazar program includes discov- 
ery observations, monitoring of known VHE blazars, 
and target-of-opportunity (ToO) observations of flar- 
ing sources. All VERITAS detections are supported 
by extensive simultaneous multiwavelength observa- 
tions in the optical, X-ray, and often GeV band, al- 
lowing detailed SED modeling to be carried out. 



Source 


z 


Type 


VERITAS references 


3C 66At 


0.444? 


IBL 


ATel 1753, [6, L3] 


1ES 0229+200 


0.140 


HBL 


[31] 


RBS0413t 


0.190 


HBL 


ATel 2272, [32] 


1ES 0414+009 


0.287 


HBL 


[31] 


1ES 0502+675+ 


0.341? 


HBL 


ATel 2301 


VERJ0521+21lt 


? 


HBL 


ATel 2260, 2309, [33] 


RXJ0648.7+1516t 


0.179 


HBL 


ATel 2486, [341 


RGB J0710+59lt 


0.125 


HBL 


ATel 1941, [351 


1ES 0806+524+ 


0.138 


HBL 


ATel 1415, [361 


Mrk421 


0.030 


HBL 


[37, 38J 


B2 1215+30 


0.130? 


IBL 


[31] 


1ES 1218+304 


0.184 


HBL 


[39, 40] 


W Comae t 


0.102 


IBL 


ATel 1422, 1565, [21, 41] 


PKS 1424+240+ 


? 


IBL 


ATel 2084, [42] 


H 1426+428 


0.129 


HBL 


[31] 


1ES 1440+1221" 


0.162 


IBL 


ATel 2786, [43] 


PG 1553+113 


0.5? 


HBL 


[44] 


Mrk501 


0.034 


HBL 


[45] 


1ES 1959+650 


0.047 


HBL 


[31] 


BL Lacertae 


0.069 


LBL 


ATel 3459, [14] 


1ES 2344+514 


0.044 


HBL 


[46] 



TABLE II: List of the 21 blazars detected at VHE with 
VERITAS, ordered by right ascension. The 10 VHE dis- 
coveries are marked with t. Unknown and uncertain red- 
shifts are indicated with a question mark (?). The blazar 
type classifications are taken from [4T| or are derived from 
VERITAS-led multiwavelength studies. 



Apart from the positive detections, VERITAS has 
observed a sample of 47 X-ray/EGRET-selcctcd tar- 
gets where no significant signal was detected. Al- 
though none of these blazars was detected individu- 
ally, stacking the results from these observations yields 
an excess with 4.1er significance. Interestingly, the 
same stacking procedure applied to 21 .Fermi-selected 
blazar targets does not show a positive excess [43| . 

VERITAS also has an extensive ToO program to 
carry out follow-up VHE observations on blazars that 
are seen in a high state at lower frequencies. VER- 
ITAS currently has an existing program using pub- 
lic Fermi-LAT data to identify high GeV gamma-ray 
states [43]. VERITAS also responds to public an- 
nouncements of blazar flares in the gamma-ray, X 
ray, and optical bands, and has developed an inter- 
nal blazar monitoring program at optical frequencies 
at the UCO/Lick Observatory. A summary of recent 
results from the ToO program can be found in [49| . 



III. RECENT BLAZAR DISCOVERIES 

The discovery of new VHE-emitting blazars is one 
of the science goals of the blazar observing program. 
Five blazar discoveries have been reported by the 
VERITAS collaboration in the last two years. 
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FIG. 4: Nightly light curve of Markarian 421 during the VERITAS 2009-2010 monitoring campaign. The source was 
detected at an average flux level of approximatively 2 Crab Nebula units over the entire season. On 2010 February 17, 
the source was observed at a flux level of ~ 8 Crab Nebula units. A zoom of the intranight light curve of the flaring 
event in 2-minute bins is shown. Taken from [381 ]. 



RBS 0413 is an HBL with a redshift of z = 0.190. It 

was a VHE candidate due to its high X-ray flux [3(. 
Deep VERITAS observations were motivated by its 
hard sjgectrum in the GeV band measured by Fermi- 
LAT 4]. RBS 0413 was observed by VERITAS for 
~26 h of quality-selected live time between September 
2008 and January 2010. An excess with statistical 
significance of 5.5er was observed. The time-averaged 
spectrum follows a power-law with photon index of 
T = 3.2 ± 0.7 and integral flux of -2% Crab Nebula 

a 

1ES 0502+675 observations were motivated by the 
flux and spectrum reported by Fermi-LAT . It was ob- 
served by VERITAS for ~30 h of quality-selected live 
time between September 2009 and January 2010. A 
gamma-ray excess corresponding to a statistical sig- 
nificance of ^ller was observed. The VHE flux is 
constant within the observed statistics, with an av- 
erage value of Crab Nebula, and a power-law 
spectrum with photon index T = 3.92 ±0.35 (43[. The 
reported redshift of 0.341 is uncertain. Following the 
VERITAS detection, optical spectroscopy was done 
to determine the redshift value, but no significant ab- 
sorption or emission lines were seen. 

VER J0521+211 was suggested as a VHE candi- 
date following a search carried out inside the VER- 
ITAS collaboration for clusters of high-energy pho- 
tons collected by Fermi-LAT. VERITAS observations 
were taken between 22 and 24 October 2009 (MJD 
55126-55128), leading to the detection of a new VHE 
gamma-ray source: VER J0521+211 (33|. An excess 
with significance of 5.5<r was obtained during the first 
230 minutes of exposure. VERITAS continued moni- 
toring the new gamma-ray source and detected a high 
flux state in 22 November 2009 (MJD 55157), when 
VER J0521+211 reached a peak flux more than three 



times as high as during the discovery observations. A 
total exposure of 14.5 hours was accumulated, leading 
to a signal with significance of 15.6<7 and power-law 
spectrum with photon index V — 3.44 ± 0.20. Follow- 
ing the VHE detection, optical spectroscopy observa- 
tions revealed a continuum-dominated spectrum that 
identifies the optical counterpart as a BL Lac-type 
blazar. No spectral absorption lines could be iden- 
tified, and therefore the redshift of VER J0521+211 
remains unknown (33j |. 

RX J0648.7+1516 is a ROSAT-detected X-ray 
source associated with a compact flat-spectrum ra- 
dio source, located 6.3° off the galactic plane. The 
source was identified as a promising VHE candidate 
by the Fermi-LAT Collaboration, and this informa- 
tion prompted VERITAS observations of the source. 
RX J0648. 7+1516 was observed with VERITAS be- 
tween 4 March and 15 April 2010 (MJD 55259-55301) 
for a total exposure of 19.3 hours of quality-selected 
live-time. A gamma-ray signal with 5.2er significance 
was detected, with a spectrum following a power law 
with photon index T = 4.4 ± 0.8 and integral flux 
of 7% that of the Crab Nebula [34J . Optical spec- 



troscopy measurements were triggered following the 
VERITAS detection. The optical spectrum showed 
Ca H+K, G-band, Na I and Mg I spectral absorption 
features, compatible with the source being at a red- 
shift of z = 0.179 [13]. These observations also provide 
a definitive optical identification of RX J0648. 7+1516 
as a blazar of the BL Lac subclass. 
1ES 1440+122 is an IBL with a redshift of z = 
0.162. This blazar was observed by VERITAS for 
~47 h between May 2008 and June 2010. An excess 
with statistical significance of 5.5cr was detected, cor- 
responding to an integral flux of ^1% Crab Nebula 
[H| with photon index T = 3.4 ± 0.7. 
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FIG. 5: Observed gamma-ray light curve of BL Lacertae 
on the night of 2011 June 28. Clear flux variability was 
seen during the 40 minutes of VERITAS observations, with 
a flux-halving time of ~ 15 min. This is the shortest vari- 
ability timescale observed so far in an LBL. 



IV. OBSERVATIONS OF FLARING 
BLAZARS 



Mrk 421 is regularly monitored by VERITAS, taking 
short snapshots at irregular intervals. On 2010 Febru- 
ary 17, VERITAS detected a very bright gamma-ray 
flare from Mrk 421 and therefore extended the cover- 
age to five hours. Figure |4] shows the Mrk 421 nightly 
light curve over the entire 2009-10 season, with a de- 
tailed view of the 2010 February 17 flaring event. Dur- 
ing the flare, the average gamma-ray flux exceeded 8 
Crab Nebula units, with a significance exceeding 260<r 
over the 4.9 h of observation time [38[ . Flux variability 
is observed on time scales shorter than 10 min. 

BL Lacertae was discovered at VHE by MAGIC dur- 
ing a flare in 2005 with a flux of ~ 3% Crab Nebula 
|50{ . The source has been repeatedly observed with 
ground-based Cherenkov telescopes since then, but it 
was not detected again until 2011. In May 2011, sev- 
eral observatories reported an increase in activity of 
the source from optical to GeV energies. VERITAS 
monitoring observations started on 2011 May 26 and 
continued on several nights through the end of June. 
On 2011 June 28, a flare was observed by VERITAS 
during a 40-minute observation under moonlight con- 
ditions. The observed gamma-ray flux exceeded 50% 
Crab Nebula during the flare, with a total significance 
above 20cr [l4[ ■ Figure [5] shows the VHE gamma-ray 
light curve during the observed flare. Clear intranight 
variability is seen, with halving time of ~ 15 min, 
which is the fastest variability ever observed in an 
LBL. 



V. DISCUSSION AND CONCLUSIONS 

The population of VHE-dctccted blazars is rapidly 
growing. The VHE blazar sample is dominated by 
HBLs, but thanks to the increased sensitivity of 
current ground-based Cherenkov telescopes it also 
includes a small number of low-frcquency-peaked 
BL Lacs and FSRQs. VERITAS has contributed 
to this progress with the discovery of 10 new VHE 
blazars in the past few years and several key observa- 
tions of known VHE sources in flaring states, always 
accompanied by simultaneous multiwavelength cover- 
age. 

The VERITAS blazar science program has several 
goals: 

• Monitoring of a sample of known nearby VHE 
blazars to study the properties of their gamma- 
ray emission: flux and spectral variability, duty 
cycle, extension of the spectrum above lOTeV, 
and measurement of the shortest variability 
timescales. 

• Extensive observations of known hard-spectrum 
VHE blazars with the main goal of character- 
izing their emission above 1 TeV to study the 
interaction of multi-TeV photons with the EBL, 
thereby constraining the EBL density and the 
magnitude of the primordial intergalactic mag- 
netic fields. 

• Monitoring of bright IBLs and LBLs to char- 
acterize their gamma-ray emission during non- 
flaring episodes and understand under which 
conditions low-frequency-peaked BL Lacs emit 
significant levels of VHE emission. 

• Coordination of broadband multiwavelength ob- 
servations simultaneous with VERITAS detec- 
tions of blazars, focusing on radio, optical, X- 
ray and GeV bands. Multiwavelength coverage 
is used to build simultaneous SEDs and perform 
correlated variability studies that can be com- 
pared with the predictions from different emis- 
sion models. 

• Discovery observations of new VHE candidates 
to enlarge the population of VHE blazars, with 
special emphasis on increasing the number of de- 
tected IBLs, LBLs and FSRQs. 

• A ToO program to follows up on reports of high 
radio/optical/X-ray/GeV states for selected tar- 
gets, in order to characterize the VHE emission 
of these objects and to study correlated flaring 
behavior between different energy bands. 

VERITAS will continue to develop a strong pro- 
gram of blazar observations in the coming years. 
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